Real Spectra in Non-Hermitian Hamiltonians Having VT Symmetry 
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The condition of self-adjointness ensures that the eigenvalues of a Hamiltonian are real and bounded 
below. Replacing this condition by the weaker condition of VT symmetry, one obtains new infinite 
classes of complex Hamiltonians whose spectra are also real and positive. These VT symmetric 
theories may be viewed as analytic continuations of conventional theories from real to complex 
phase space. This paper describes the unusual classical and quantum properties of these theories. 
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Several years ago, D. Bessis conjectured on the basis of 
numerical studies that the spectrum of the Hamiltonian 
H = + + ix^ is real and positive . To date there is 
no rigorous proof of this conjecture. We claim that the re- 
ality of the spectrum of H is due to VT symmetry. Note 
that H is invariant neither under parity V, whose effect 
is to make spatial reflections, p —p and x — s- —x, nor 
under time reversal T, which replaces p —p, x x, 
and i — > —i. However, VT symmetry is crucial. For ex- 



ample, the Hamiltonian p 



IX 



ix has VT symmetry 



and our numerical studies indicate that its entire spec- 
trum is positive definite; the Hamiltonian p^ -f ix^ -I- a; is 
not PT-symmetric, and the entire spectrum is complex. 

The connection between VT symmetry and positivity 
of spectra is simply illustrated by the harmonic oscillator 
H = p'^+x'^, whose energy levels are En = 2n+l. Adding 
ix to H does not break VT symmetry, and the spectrum 
remains positive definite: En = 2n + j. Adding —a; also 
does not break VT symmetry if we define V as reflection 
about X = ■^j X ^ 1 — X, and again the spectrum remains 
positive definite: En = 2n + |. By contrast, adding 
ix — X does break VT symmetry, and the spectrum is 
now complex: En = 2n + 1 + ■^i. 

The Hamiltonian studied by Bessis is just one example 
of a huge and remarkable class of non-Hermitian Hamil- 
tonians whose energy levels are real and positive. The 
purpose of this Letter is to understand the fundamen- 
tal properties of such a theory by examining the class of 
quantum-mechanical Hamiltonians 



H 
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p + m X 



[ixy^ {N real). 



(1) 



As a function of N and mass we find various phases 
with transition points at which entirely real spectra begin 
to develop complex eigenvalues. 

There are many applications of non-Hermitian VT- 
invariant Hamiltonians in physics. Hamiltonians ren- 
dered non-Hermitian by an imaginary external field have 
been introduced recently to study delocalization transi- 
tions in condensed matter systems such as vortex flux- 
line depinning in type-II superconductors or even to 
study population biology H]. Here, initially real eigenval- 



ues bifurcate into the complex plane due to the increas- 
ing external field, indicating the unbinding of vortices 
or the growth of populations. We believe that one can 
also induce dynamic delocalization by tuning a physical 
parameter (here N) in a self-interacting theory. 

Furthermore, it was found that quantum field the- 
ories analogous to the quantum-mechanical theory in 
Eq. (|l|) have astonishing properties. The Lagrangian 
L — (y (f))^ +m?<lP — g{i(f))^ {N real) possesses PT invari- 
ance, the fundamental symmetry of local self-interacting 
scalar quantum field theory [Q. Although this theory 
has a non-Hermitian Hamiltonian, the spectrum of the 
theory appears to be positive definite. Also, L is explic- 
itly not parity invariant, so the expectation value of the 
field {(j)) is nonzero, even when AT = 4 |]. Thus, one 
can calculate directly (using the Schwinger-Dyson equa- 
tions, for example |^) the (real positive) Higgs mass in 
a renormalizable theory such as — or igij)'^ in which 
symmetry breaking occurs naturally (without introduc- 
ing a symmetry-breaking parameter) . 

Replacing conventional gcj)'^ or gcj)^ theories by —gcf)^ 
or ig(jfi theories has the effect of reversing signs in the 
beta function. Thus, theories that are not asymptotically 
free become asymptotically free and theories that lack 
stable critical points develop such points. For example, 
T^T-symmetric masslcss electrodynamics has a nontrivial 
stable critical value of the fine-structure constant a [0 . 

Supersymmetric non-Hermitian, T^T-invariant La- 
grangians have been examined ||]. It is found that the 
breaking of parity symmetry does not induce a break- 
ing of the apparently robust global supersymmetry. The 
strong-coupling limit of non-Hermitian PT-symmetric 
quantum field theories has been investigated |^ ; the cor- 
related limit in which the bare coupling constants g and 
— both tend to infinity with the renormalized mass 
M held fixed and finite, is dominated by solitons. (In 
parity-symmetric theories the corresponding limit, called 
the Ising limit, is dominated by instantons.) 

To elucidate the origin of such novel features we ex- 
amine the elementary Hamiltonian (|^) using extensive 
numerical and asymptotic studies. As shown in Fig. 
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when m = the spectrum of H exhibits three distinct 
behaviors as a function of TV. When N > 2, the spectrum 
is infinite, discrete, and entirely real and positive. (This 
region includes the case N = A for which H = — x'^; 
the spectrum of this Hamiltonian is positive and discrete 
and (a:) ^ in the ground state because H breaks parity 
symmetry!) At the lower bound = 2 of this region 
lies the harmonic oscillator. A phase transition occurs 
at A = 2; when 1 < A^ < 2, there are only a finite 
number of real positive eigenvalues and an infinite num- 
ber of complex conjugate pairs of eigenvalues. In this 
region VT symmetry is spontaneously broken . As A^ 
decreases from 2 to 1, adjacent energy levels merge into 
complex conjugate pairs beginning at the high end of the 
spectrum; ultimately, the only remaining real eigenvalue 
is the ground-state energy, which diverges as A^ ^ 1+ 
p| . When A^ < 1, there are no real eigenvalues. The 
massive case m ^ is even more elaborate; there is a 
phase transition at A^ = 1 in addition to that at A = 2. 
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FIG. 1. Energy levels of the Hamiltonian H = ~ {i^)^ 
as a function of the parameter N . There are three regions: 
When N > 2 the spectrum is real and positive. The lower 
bound of this region, N = 2, corresponds to the harmonic 
oscillator, whose energy levels are E„ = 2n + 1. When 
1 < A < 2, there are a finite number of real positive eigen- 
values and an infinite number of complex conjugate pairs of 
eigenvalues. As N decreases from 2 to 1, the number of real 
eigenvalues decreases; when N < 1.42207, the only real eigen- 
value is the ground-state energy. As N approaches 1+, the 
ground-state energy diverges. For A < 1 there are no real 
eigenvalues. 

The Schrodinger eigenvalue differential equation corre- 
sponding to the Hamiltonian (|l|) with to = is 

-il;"{x) - {ix)^il;{x) = Eil;{x). (2) 

Ordinarily, the boundary conditions that give quantized 



energy levels E are that ip{x) ^ as |x| — > oo on the 
real axis; this condition suffices when 1 < A^ < 4. How- 
ever, for arbitrary real A^ we must continue the eigenvalue 
problem for into the complex- a: plane. Thus, we re- 
place the real-x axis by a contour in the complex plane 
along which the differential equation holds and we im- 
pose the boundary conditions that lead to quantization 
at the endpoints of this contour. (Eigenvalue problems 
on complex contours are discussed in Ref. Q.) 




FIG. 2. Wedges in the complex- a: plane containing the 
contour on which the eigenvalue problem for the differential 
equation (^) for = 4.2 is posed. In these wedges il){x) 
vanishes exponentially as — > oo. The wedges are bounded 
by Stokes lines of the differential equation. The center of the 
wedge, where V'(^) vanishes most rapidly, is an anti-Stokes 
line. 

The regions in the cut complex- a; plane in which ip{x) 
vanishes exponentially as |x| oo are wedges (see Fig.||); 
these wedges are bounded by the Stokes lines of the dif- 
ferential equation jl3|. The center of the wedge, where 
ijj{x) vanishes most rapidly, is called an anti-Stokes line. 

There are many wedges in which ip{x) — > as |x| — s- cx). 
Thus, there are many eigenvalue problems associated 
with a given differential equation jl2|. However, we 
choose to continue the eigenvalue equation (|^) away from 
the conventional harmonic oscillator problem at A^ = 2. 
The wave function for N — 2 vanishes in wedges of angu- 
lar opening ^tt centered about the negative- and positive- 
real X axes. For arbitrary A^ the anti-Stokes lines at the 
centers of the left and right wedges lie at the angles 



A^-2 vr 



71oft = -TT - 



A^ + 2 2 



and bright = - 



A^-2 TT 



A^ + 2 2 



(3) 



The opening angle of these wedges is A = 2tt/{N + 2). 
The differential equation (|^) may be integrated on any 
path in the complex- a; plane so long as the ends of the 
path approach complex infinity inside the left wedge and 



2 



the right wedge . Note that these wedges contain the 
real- a; axis when 1 < iV < 4. 

As N increases from 2, the left and right wedges rotate 
downward into the complex- a; plane and become thinner. 
At A'^ = oo, the differential equation contour runs up 
and down the negative imaginary axis and thus there is 
no eigenvalue problem at all. Indeed, Fig. |l| shows that 
the eigenvalues all diverge as A^ ^ cxd. As A^ decreases 
below 2 the wedges become wider and rotate into the 
upper-half x plane. At A = 1 the angular opening of 
the wedges is |7r and the wedges are centered at |7r and 
ivr. Thus, the wedges become contiguous at the positive- 
imaginary X axis, and the differential equation contour 
can be pushed off to infinity. Consequently, there is no 
eigenvalue problem when A^ = 1 and, as we would expect, 
the ground-state energy diverges as N ^1+ (see Fig.l). 

To ensure the numerical accuracy of the eigenvalues 
in Fig. |l|, we have solved the differential equation (||) 
using two independent procedures. The most accurate 
and direct method is to convert the complex differential 
equation to a system of coupled, real, second-order equa- 
tions which we solve using the Runge-Kutta method; the 
convergence is most rapid when we integrate along anti- 
Stokes lines. We then patch the two solutions together at 
the origin. We have verified those results by diagonaliz- 
ing a truncated matrix representation of the Hamiltonian 
in Eq. in harmonic oscillator basis functions. 

Semiclassical analysis: Several features of Fig. ^ can be 
verified analytically. When N > 2, WKB gives an excel- 
lent approximation to the spectrum. The novelty of this 
WKB calculation is that it must be performed in the 
complex plane. The turning points x± are those roots of 
E + (ix)^ — that analytically continue off the real axis 
as A^ moves away from A = 2 (the harmonic oscillator): 

X_ = ^l/JVg.^(3/2-l/JV)^ ^ ^l/N^-^.il/2-l/N)^ (4) 

These turning points lie in the lower-half (upper-half) x 
plane in Fig. || when A^ > 2 (A^ < 2). 

The leading-order WKB phase-integral quantization 
condition is (n + l/2)7r = Ji^^ dx y/WTJix)^ . It IS cru- 
cial that this integral follow a path along which the in- 
tegral is real. When N > 2, this path lies entirely in the 
lower-half x plane and when N = 2 the path lies on the 
real axis. But, when N < 2 the path is in the upper- half 
X plane; it crosses the cut on the positive-imaginary 

TABLE I. Comparison of the exact eigenvalues (obtained 
with Runge-Kutta) and the WKB result in (H). 



N 


n 


exact 


-EwKB 


TV 


n 


-E^cxact 


-EwKB 


3 





f.l562 


1.0942 


4 





1.4771 


1.3765 




1 


4.1092 


4.0894 




1 


6.0033 


5.9558 




2 


7.5621 


7.5489 




2 


11.8023 


11.7689 




3 


11.3143 


11.3042 




3 


18.4590 


18.4321 




4 


15.2916 


15.2832 











axis and thus is not a continuous path joining the turning 
points. Hence, WKB fails when N < 2. 

When A^ > 2, we deform the phase-integral contour so 
that it follows the rays from x^ to and from to 2:+: 
{n + 1/2)71 ^ 2sm{TT/N)E^/^+^/^ dsVT 
then solve for E„ : 



We 



E„ 



r(3/2 + l/A^)0F(n + 1/2) 



sin(7r/Ar)r(l + 1/N) 



2N 
N + 2 



in 



3). (5) 



We perform a higher-order WKB calculation by replacing 
the phase integral by a closed contour that encircles the 
path in Fig. | (see Ref. |l|,0). See Table I. 

It is interesting that the spectrum of the potential 
is like that of the —{ix)^ potential. The leading-order 
WKB quantization condition (accurate for A^ > 0) is 
like Eq. (||) except that sin(7r/A^) is absent. However, 
as A^ — > 00, the spectrum of approaches that of 
the square- well potential = {n + l)^7r^/4], while the 



energies of the 



potential diverge (see Fig. 1). 



Asymptotic study of the ground-state energy near N = 1: 
When A = 1, the differential equation (|^) can be solved 
exactly in terms of Airy functions. The anti-Stokes lines 
at A^ = 1 lie at 30° and at 150°. We find the solution 
that vanishes exponentially along each of these rays and 
then rotate back to the real- a; axis to obtain 

V'left, right (x) = Ci,2 Ai(Ta:e±''^/6 + i?e±2-/3). (6) 

We must patch these solutions together at a; = accord- 
ing to the patching condition ^l'0(2^)P|^_o ~ ^' 
real E, the Wronskian identity for the Airy function is 



dx 



|Ai(xe-*'^/6_,_^g-2»^/3)|5 



x=0 



1 

'2^ 



(7) 



instead of 0. Hence, there is no real eigenvalue. 

Next, we perform an asymptotic analysis for A = 1 + e, 
—il}"{x) — {ixy~^'^tp{x) = E^/j{x), and take V'(a;) =1/0(2:) + 
eyi{x) -f O(e^) as e ^ 0+. We assume that — > 00 as 
e 0-1- , let C2 = 1 in Eq. (^, and obtain 



yo 



(0) = Ai(i;e-2W3) _ e'-/eE-^/^ei^"'"/{2V^). (8) 



We set 2/1 (x) = Q{x)yQ{x) in the inhomogeneous equation 
—y'{{x) — ixyi(x) — Eyi{x) = ix\n.{ix)yQ{x) and get 



?'(0) 



vm Jo 



dx X \n{ix)yQ{x). 



(9) 



Choosing Q{0) = 0, we find that the patching con- 
dition at a; = gives 1 = 27re |?/o(0)|^ [Q'(0) + Q*'(0)], 
where we have used the zeroth-order result in Eq. ([^). 
Using Eqs. (^ and (||) this equation becomes 



1 = 



Re 
E 



vm Jo 



dxx \n{ix)yQ{x) 



(10) 
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Since j/o(a;) decays rapidly as x increases, the integral in 
Eq. (|l^) is dominated by contributions near 0. Asymp- 
totic analysis of this integral gives an implicit equation 
for E' as a function of e (see Table II): 

1 - ee^^'^'s-3/2[V3in(2V;B) + tt - (1 - -f)V3]/8. (11) 

Behavior near N — 2: The most interesting aspect of 
Fig. |l] is the transition that occurs at = 2. To de- 
scribe quantitatively the merging of eigenvalues that be- 
gins when < 2 we let iV = 2 — e and study the asymp- 
totic behavior as e — > 0-f . (A Hermitian perturbation 
causes adjacent energy levels to repel, but in this case 
the non-Hermitian perturbation of the harmonic oscil- 
lator {ix)^~'^ ~ — ex^[ln(|x| -|- ^Z7rsgn(x)] causes the 
levels to merge.) A complete description of this asymp- 
totic study is given elsewhere [|lO| . 

The onset of eigenvalue merging is a phase transition 
that occurs even at the classical level. Consider the clas- 
sical equations of motion for a particle of energy E sub- 
ject to the complex forces described by the Hamiltonian 
(|l|). For TO = the trajectory x{t) of the particle obeys 
±dx[E + (ia;)^]"i/2 ^ 2dt. While E and dt are real, 
x(t) is a path in the complex plane in Fig. |^; this path 
terminates at the classical turning points x± in (^ . 

When N > 2, the trajectory is an arc joining x± in the 
lower complex plane. The motion is periodic; we have a 
complex pendulum whose (real) period T is 



T 



2- Af 

2E 2iv cos 



(A^- 2)7r' 
27V 



r(i + i/A^)V^ 

T{l/2 + l/N) ■ 



(12) 



At A^ = 2 there is a global change. For A^ < 2 a 
path starting at one turning point, say x-^-, moves to- 
ward but misses the turning point X- . This path spirals 
outward crossing from sheet to sheet on the Riemann 
surface, and eventually veers off to infinity asymptotic 
to the angle 2^v^" Hence, the period abruptly becomes 
infinite. The total angular rotation of the spiral is finite 
for all A^ 7^ 2 and as A^ ^ 2+, but becomes infinite as 
N 2~ . The path passes many turning points as it spi- 
rals anticlockwise from x+. [The nth turning point lies 
at the angle ^" (^+ corresponds to ti = 0).] As A^ 

approaches 2 from below, when the classical trajectory 
passes a new 

TABLE II. Comparison of the exact ground- state energy 
E near A = 1 and the asymptotic results in Eq. (|ll|). The 
explicit dependence of _E on e is roughly E oc { — \ne)^^^ . 



N ~ 1 



Eq. (11) 



0.1 

0.01 

0.001 

0.0001 

0.00001 

0.000001 

0.0000001 



1.6837 
2.6797 
3.4947 
4.1753 
4.7798 
5.3383 
5.8943 



2.0955 
2.9624 
3.6723 
4.3013 
4.8776 
5.4158 
5.9244 



turning point, there corresponds an additional merging 
of the quantum energy levels as shown in Fig. |^). This 
correspondence becomes exact in the limit N 2~ and 
is a manifestation of Ehrenfest's theorem. 
Massive case: The m ^ analog of Fig. |l| exhibits a new 
transition at A^ = 1 (see Fig. ||). As A^ approaches 1 
from above, the energy levels reemerge from the complex 
plane in pairs and at A^ = 1 the spectrum is again entirely 
real and positive. Below A^ = 1 the energies once again 
disappear in pairs, now including the ground state. As 
A^ — > the infinite real spectrum reappears again. The 
massive case is discussed further in Ref. [|lO|. 
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FIG. 3. The m 7^ analog of Fig. 1. Note that transitions 
occur at A = 2 and N = 1. 
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